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Abstract 


The earliest confirmed interstellar object, 'Oumuamua, was discovered in the solar system by Pan-STARRS in 
2017, allowing for a calibration of the abundance of interstellar objects of its size ~100 m. This was followed by 
the discovery of Borisov, which allowed for a similar calibration of its size ~0.4—1 km. One would expect a much 
higher abundance of significantly smaller interstellar objects, with some of them colliding with Earth frequently 
enough to be noticeable. Based on the CNEOS catalog of bolide events, we identify the ~0.45 m meteor detected 
at 2014 January 8 17:05:34 UTC as originating from an unbound hyperbolic orbit. The U.S. Department of 
Defense has released an official letter stating that “the velocity estimate reported to NASA is sufficiently accurate 
to indicate an interstellar trajectory,” which we rely on here as confirmation of the object’s interstellar trajectory. 
Based on the data provided by CNEOS, we infer that the meteor had an asymptotic speed of 
Voo ~ 42.1 3: 5.5 kms! outside of the solar system. Note that væ here refers to the velocity of the meteor 
outside the solar system, not the velocity of the meteor outside the atmosphere. Its origin is approximately toward 
R.A. 49°4 + 4°1 and decl. 11?2 + 1°8, implying that its initial velocity vector was 58 + 6 km s™! away from the 


velocity of the local standard of rest (LSR). 


Unified Astronomy Thesaurus concepts: Meteors (1041); Interstellar objects (52) 


1. Introduction 


*Oumuamua was the first interstellar object detected in the 
solar system by Pan-STARRS (Meech et al. 2017; Micheli et al. 
2018). Several follow-up studies of 'Oumuamua were conducted 
to better understand its origin and composition (Bannister et al. 
2017; Bolin et al. 2017; Gaidos et al. 2017; Jewitt et al. 2017; 
Mamajek 2017; Ye et al. 2017; Bialy & Loeb 2018; 
Fitzsimmons et al. 2018; Hoang et al. 2018; Trilling et al. 
2018; Seligman et al. 2019; Siraj & Loeb 2019a, 2019b). Its size 
was estimated to be 20-200m, based on Spitzer Space 
Telescope constraints on its infrared emission given its 
temperature (Trilling et al. 2018). The discovery of 'Oumuamua 
was followed by that of the second interstellar object, Borisov, in 
2019 (Guzik et al. 2020). The size of Borisov's nucleus was 
estimated to be 0.4—1 km (Jewitt et al. 2020; Siraj & Loeb 2021). 

Forbes & Loeb (2019) predicted that spectroscopy of 
*Oumuamua-like objects grazing the Sun could reveal their 
chemical compositions. Since there should be a higher 
abundance of interstellar objects smaller than ‘Oumuamua, 
we could observe small interstellar objects impacting the 
Earth's atmosphere, appearing as interstellar meteors (Opik 
1950). Spectroscopy of the gaseous debris from such objects as 
they burn up in the Earth's atmosphere could reveal their 
composition. There is significant evidence for previous 
detections of dust-sized interstellar meteors (Baggaley et al. 
1993; Taylor et al. 1996; Mathews et al. 1999; Baggaley 2000; 
Meisel et al. 2002; Weryk & Brown 2004; Afanasiev et al. 
2007; Froncisz et al. 2020), as well as interstellar meteor 
nondetections (Hajdukova 1994; Musci et al. 2012; Engelhardt 
et al. 2017; Hajdukova et al. 2020) but to date no definitive 
evidence of any meter-scale interstellar meteors. 


Original content from this work may be used under the terms 

BY of the Creative Commons Attribution 4.0 licence. Any further 
distribution of this work must maintain attribution to the author(s) and the title 
of the work, journal citation and DOI. 


The CNEOS catalog includes the geocentric velocity 
components and geographic coordinates for bolides detected 
by U.S. government sensors. In this article, we identify a 
meteor from the CNEOS catalog that is likely of interstellar 
origin. Furthermore, we explore the size distribution of 
interstellar meteors and motivate the investigation of interstellar 
meteors. We present a strategy for conducting spectroscopy 
and obtaining physical samples of interstellar meteors. 


2. Methods 


We analyzed the bolide events in the CNEOS catalog, and 
found that the meteor detected at 2014 January 8 17:05:34 UTC 
had an unusually high heliocentric velocity at impact. 
Accounting for the motion of the Earth relative to the Sun 
and the motion of the meteor relative to the Earth, we found 
that the meteor had a heliocentric velocity of ~60kms_' at 
impact, which implies that the object was unbound. To uncover 
the kinematic history of this meteor, we integrated its motion 
from impact backward in time. 

The Python code created for this work used the open-source 
N-body integrator software REBOUND’ to trace the motion of 
the meteor under the gravitational influence of the solar system 
(Rein & Liu 2012). 

We initialize the simulation with the Sun, the eight planets, 
and the meteor, with geocentric velocity vector (VXobs» VYobss 
VZops) = ( —3.4, —43.5, — 10.3) kms’, located at 1°3 S 147°6 
E, at an altitude of 18.7 km, at the time of impact, t; = 2014 
January 8 17:05:34 UTC, as reported in the CNEOS catalog. 
We then use the IASI5 adaptive time-step integrator to trace 
the meteor’s motion back in time (Rein & Spiegel 2014). This 


l https: //cneos.jpl.nasa.gov /fireballs / 

? The fastest meteor in the CNEOS catalog obtains its high speed from a 
head-on orbit relative to the Earth and its extrapolated orbit is found to be 
bound to the Sun. The meteor we focus on is the second fastest. The orbit of the 
third fastest meteor in the catalog is possibly bound within uncertainties. 

* https: / [rebound.readthedocs.io /en/latest / 
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does not account for air drag, which would lead to an even 
higher impact speed, and therefore heliocentric speed, given the 
encounter geometry. The slowdown due to air drag is estimated 
in Siraj & Loeb (2022). 


3. Results 
3.1. Trajectory 


There are no substantial gravitational interactions between 
the meteor and any planet other than Earth for any trajectory 
within the reported errors. Based on the impact speed reported 
by CNEOS, vobs = 44.8 km s !, we find that the meteor was 
unbound with an asymptotic speed of v4, 42.1 kms"! 
outside of the solar system. In order for the object to be 
bound, the observed speed of v,,, = 44.8 km s ! would have to 
be off by more than 45%, or 20 kms |, or assuming a correct 
speed, a radiant off by more than 60? (Zuluaga 2019). 

The typical velocity uncertainty for meter-scale impactors in 
the CNEOS catalog was estimated by Brown et al. (2016) and 
Granvik & Brown (2018) to be less than 1 km s |, but 2 of 10 
events analyzed by Devillepoix et al. (2019) have discrepancies 
with alternate measurements of up to 2896 in speed. Drolshagen 
et al. (2020) noted that the CNEOS catalog includes two 
objects that were detected in space before impact, 2008 TC; 
and 2018 LA, with respective measurement errors of 0.3 and 
0.52 kms ' on the geocentric impact speed. These papers 
underlined the importance of uncertainties and motivated us to 
retrieve the measurement uncertainties for the 2014 January 8 
meteor specifically. 

Recently, the U.S. Department of Defense, which houses the 
classified data pertaining to the uncertainties involved in the 
CNEOS 2014-01-08 detection, released a public statement 
dated 2022 March 1 and addressed to the NASA Science 
Mission Directorate, referencing this discovery preprint and 
mentioning our analysis that the meteor originated from an 
unbound hyperbolic orbit with high confidence (Shaw 2022). 
The letter then states “Dr. Joel Mozer, the Chief Scientist of 
Space Operations Command, reviewed analysis of additional 
data available to the Department of Defense related to this 
finding. Dr. Mozer confirmed that the velocity estimate 
reported to NASA is sufficiently accurate to indicate an 
interstellar trajectory" (Shaw 2022). We rely on the official 
letter released by the U.S. Department of Defense for the 
confirmation of the interstellar origin of the 2014 January 8 
meteor. Given the confidence level expressed in the official 
letter from the U.S. Department of Defense, we calculate and 
quote here the corresponding Gaussian symmetric errors on 
each parameter related to the object's orbital trajectory. 

We find that the heliocentric orbital elements of the meteor at 
time of impact are as follows: semimajor axis, a= —0.47 + 
0.15 au, eccentricity, e = 2.4 + 0.3, inclination i= 10 2°, 
longitude of the ascending node, Q = 108 + 1°, argument of 
periapsis, w =58 + 2°, and true anomaly, f= —58 + 2?. The 
trajectory is shown in Figure 1. The origin is toward R.A. 
49.4 + 4.1? and decl. 11.2 + 1.8?. The heliocentric incoming 
velocity at infinity of the meteor in right-handed Galactic 
coordinates is v4,(U, V, W) = (32.7 € 5.8, —4.5 + 1.5, 26.1 + 
2.0) km s~}, which is 58 + 6 km s~! away from the velocity of 
the local standard of rest (LSR), (U, V, W)zsr — (—11.1, 
—12.2, —7.3) km s7! (Schonrich et al. 2010). 
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Figure 1. Trajectory of the 2014 January 8 meteor (red), shown intersecting 
with that of Earth (blue) at the time of impact, t; = 2014 January 8 17:05:34. 


3.2. Size Distribution 


Given the impact speed of the meteor, ~44.8 km s^, and the 
total impact energy, 4.6 x 10'8 ergs, the meteor mass was 
approximately 4.6 x 10° g. Adopting bulk density values of 1.7 
and 0.9 g cm ? for Type II and Type IIa objects respectively, 
yields a radius of 0.4-0.5m for a spherical geometry 
(Ceplecha 1988; Palotai et al. 2019). Adopting a density 
similar to an iron meteorite, ~8 g cm ^, as suggested by Siraj 
& Loeb (2022) based on the light curve for CNEOS 2014-01- 
08, yields a diameter of ~0.5 m. 

The CNEOS catalog includes bolide events at a relatively 
high frequency for the past decade, so we approximate the 
yearly detection rate of interstellar meteors to be at 
least~ 0.1 yr. This estimate agrees with that of Hajdukova 
et al. (2019) within the error budget for the size of CNEOS 
2014-01-08. We estimate the number density of similarly sized 
interstellar objects by dividing the yearly detection rate by the 
product of the impact speed of the meteor and the cross 
sectional area of the Earth, finding the approximate number 
density of interstellar objects with a size of order R~ 0.45 m 
and a speed v ~ 60kms ‘kms | relative to the LSR to be 


» 0.1 yr! 
(13 AU/yr)(5.7 x 10-?AU?) 


Given 9546 Poisson uncertainties (Gehrels 1986), the 
inferred local number density for interstellar objects of this 
size is n = 10515 au~’. This figure necessitates 6 x 102215 
similarly sized objects, or 0.2—20 Earth masses of material, to 
be ejected per local star. This is at tension with the fact that a 
minimum-mass solar nebula is expected to have about an Earth 
mass of total planetesimal material interior to the radius where 
the orbital speed is ~60 km s ! (Desch 2007), with similar 
values for other planetary systems (Kuchner 2004). Our 
inferred abundance for interstellar meteors should be viewed 


—1065AU3, (1) 


^ Gravitational focusing by the Earth is negligible since the meteor speed 
exceeds considerably the escape speed from the Earth. The density 
enhancement due to gravitational focusing by the Sun is well below the 
uncertainty in the estimated value of n, so that our inferred range of local values 
also corresponds to the density outside of the solar system. 
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Figure 2. CNEOS 2014-01-08-like objects, ‘Oumuamua-like objects, and 
Borisov-like objects in size-abundance parameter space, expressed as number 
per star per differential unit of log size. The vertical error bars correspond to the 
95% Poisson uncertainties, while the horizontal error bars for ‘Oumuamua and 
Borisov correspond to the discrete ranges discussed in the text. 


as a lower limit since the CNEOS data might have a bias 
against detection of faster meteors (Brown et al. 2016). 

CNEOS 2014-01-08, ‘Oumuamua, and Borisov together serve 
as important calibration points for the size distribution of 
interstellar objects. Figure 2 illustrates the three classes of 
interstellar objects considered here in size-abundance parameter 
space. We use a Monte Carlo simulation to characterize the slope 
of this size distribution. For each run of the simulation, we select a 
size and a population abundance for the three interstellar objects. 
The sizes for the CNEOS 2014-01-08, ‘Oumuamua, and Borisov 
are drawn from the triangular distributions spanning 0.3-0.7 m? 
20-200 m (Trilling et al. 2018), and 0.4-1 km (Jewitt et al. 
2020; Siraj & Loeb 2021), respectively, with peaks at 0.45 m, 
100m, and 0.7km. The abundances are drawn from the 
Poisson distributions for sample size of one with central values 
of 10°, 10^! (Do et al. 2018; Levine & Laughlin 2021; Siraj & 
Loeb 2021), and 9 x 10 ? au, respectively. Given the three 
size-abundance tuples, each run of the simulation fits a linear 
least square to the log-log transformed values, resulting in a 
power-law size distribution fit with normalization k and slope 
q, corresponding to a size distribution with the form 
N(>R)~ kR' ^4. We repeat this routine 10° times, thereby 
constructing a distribution of q corresponding to the uncertain- 
ties in CNEOS 2014-01-08-like, 'Oumuamua-like, and 
Borisov-like populations. Our model applies if 'Oumuamua 
was not primarily composed of hydrogen or helium. 

We find that, over the size range from CNEOS 2014-01-08 
(meter-scale) to Borisov (kilometer-scale), the best-fit slope is 
q ~ 3.6 + 0.5, where the quoted uncertainty corresponds to 2o. 
The result is in good agreement with q — 3.5 the analytical 
collisional model by Dohnanyi (1969), but is not inconsistent 
with the q — 4 scale-free power-law model (Siraj & Loeb 2021), 
which contains equal mass per logarithmic bin. The power-law 
extrapolation may not hold for all bolide radii down to dust 
particles. 


4. Detection Strategy 
The cores of meteoroids with radii larger than ~5 cm can reach 
the ground in the form of meteorites (Kruger & Grun 2014). 


? Encapsulating a factor of 3 uncertainty in density in both directions, as well 
as +20% uncertainty in the velocity, to account for measurement error. 
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Additionally, meteoroids on smaller size scales could be 
accelerated from the Poynting-Robertson effect (Wyatt & 
Whipple 1950) and could have potential origins in the interstellar 
medium. Hence, interstellar meteors with size 25 cm are optimal 
for a systematic study of physical extrasolar material (in addition 
to the spectroscopy of the hot gases as the meteor burns up). We 
consider their potential detection here even though they represent 
a size scale from CNEOS 2014-01-08. Since we expect interstellar 
meteors with size 25 cm to strike the Earth at least a few times per 
year based on extrapolating the size distributions considered here 
(q — 3.6 + 0.5) to objects ~5 times smaller in size than CNEOS 
2014-01-08, a network of all-sky camera systems monitoring the 
sky above all land on Earth could detect an interstellar meteor 
with size 75 cm every few years, given that the land fraction of 
Earth's surface is ~30%. Such detections can be made with 
science-grade video cameras, such as those used in AMOS, 
CAMO, and CAMS? (Jenniskens et al. 2011; Weryk et al. 2013; 
Tóth et al. 2015; Jenniskens et al. 2016). These systems could 
employ the new trajectory estimation software (Vida et al. 
2020a, 2020b) that seems to be promising. 

A conservative estimate for the total area of ~70 km altitude 
atmosphere visible from a system of two all-sky camera 
systems separated by 100 km is 5 x 10? km, so to cover all 
land on Earth would require ~300 systems, or ~600 total all- 
sky camera systems, similar to CAMO but with an all-sky field 
of view, like AMOS and CAMS. 

We therefore advocate for a network of all-sky camera 
systems to conduct real-time remote spectroscopy of the hot 
gases as r 2 5 x 10 ? m interstellar meteors burn up, and to 
precisely determine their trajectories for the immediate retrieval 
of interstellar meteorite samples. 

By extrapolating the trajectory of each meteor backward in 
time and analyzing the relative abundances of each meteor's 
chemical isotopes, one can match meteors to their parent stars 
and reveal insights into planetary system formation. 
R-processed elements, such as Eu, can be detected in the 
atmospheres of stars (Frebel et al. 2016), so their abundances in 
meteor spectra can serve as important links to parent stars. 


5. Discussion 


If interstellar objects typically have ejection speeds near zero 
from their parent stars (because they are tidally removed by 
passing stars or the Milky Way from the outskirts of their 
planetary system), we would expect interstellar objects to 
reflect the dynamics of local stars. As a result, most interstellar 
objects should originate from within the local stellar velocity 
dispersion away from the LSR (Siraj & Loeb 2020). The large 
deviation for this object from the LSR of ~58+6km s ! 
suggests that it perhaps originated in the thick disk, which has 
velocity dispersion components of (oy, oy, ow) = (50, 50, 50) 
kms relative to the LSR (Bland-Hawthorn & Gerhard 2016). 
However, the ratio of local thick disk stars to thin disk stars is 
0.04, making this a minority population. Alternatively, for a 
parent planetary system with a more typical velocity relative to 
the LSR, the object could have originated in the deep interior, 
where the orbital speeds of objects are of the necessary 
magnitude. Either way, the meteor seems to have an unusual 
origin. There is also a possibility that the parameters are 
erroneously determined. 


© http: / [ cams.seti.org / 
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Siraj & Loeb (2022) analyzed the recently released light 
curve for CNEOS 2014-01-08,’ finding that the measured 
velocity and three observed flares down to an altitude of 
18.7 km imply ambient ram pressure in the range of 113-194 
MPa when the meteor disintegrated. The required yield strength 
is 220 times higher than stony meteorites and Z2 times larger 
than iron meteorites. Additionally, the implied slowdown in the 
atmosphere suggests a significantly higher initial speed, 
strengthening the case for an interstellar origin of this meteor 
and making it an outlier relative to the velocity dispersion of 
local stars. Using the same data, Pefia-Asensio et al. (2022) 
corroborated our finding of a high material strength, quoting an 
even higher value of 222.8 + 67.6 MPa. However, the reality of 
the combination of the fireball's reported altitude and speed 
remains questionable. 

The discovery of additional interstellar meteors will serve as 
an important calibration for population-wide parameters of 
interstellar objects, including their abundance and origin. 

We estimate the impact rate of similarly sized objects with 
the Earth, given 95% Poisson distribution confidence intervals, 
to be at least i100 yr Future meteor surveys could flag 
incoming objects with excess heliocentric velocities for follow- 
up preimpact observations. Spectroscopy of gaseous debris 
from these objects as they burn up in the Earth’s atmosphere 
would reveal their composition. Given that some isotope ratios 
are expected to be markedly different for objects of interstellar 
origin compared to the solar system, one could validate an 
interstellar origin (Lingam & Loeb 2018; Forbes & Loeb 2019). 
Precision tracking with the upcoming Vera Rubin Observa- 
tory’s Legacy Survey of Space and Time (LSST®) could 
determine the trajectory of meteors of interstellar origin to their 
parent systems in the Gaia catalog.’ Our discovery also implies 
that at least 4.5 x 10815 similarly sized interstellar bolide 
events have occurred over Earth’s lifetime. Potentially, 
interstellar meteors could deliver life from another planetary 
system and mediate panspermia (Ginsburg et al. 2018). 
Interestingly, the high speed for the meteor discussed here 
implies a likely origin in the habitable zone of the abundant 
population of dwarf stars, indicating that similar objects could 
carry life from their parent planetary systems. Note that this 
discussion of panspermia and the habitable zone in connection 
with CNEOS 2014-01-08 is only speculation. 

We identified one other interstellar meteor candidate from 
the entire catalog, a bolide that struck (40.5N, 18.0W) at 
04:16:37 UTC on 2017 March 9, and appeared to have a 
heliocentric speed of ~5Okms~'. We advocate for the 
declassification of the uncertainties associated with the USG 
measurements for CNEOS 2017-03-09 so its interstellar nature 
can be confirmed. The apparent lack of a set of velocities that 
disperse continuously beyond the parabolic limit is consistent 
with the notion that the CNEOS catalog is of high quality for 
the determination of interstellar meteors. Pefía-Asensio et al. 
(2022) display the dynamical properties of all of the fireballs 
reported in the CNEOS catalog in Figure 2, and enumerate the 
velocities for the apparent hyperbolic orbits in Table 8. Three 
of the five objects with apparent hyperbolic orbits have 
heliocentric speeds within <5% of the hyperbolic threshold, 
and the remaining are CNEOS 2014-01-08 and the other bolide 


Í https: / /cneos.jpl.nasa.gov /fireballs /Ic /bolide.2014.008.170534.pdf 
8 A 
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? https: //gea.esac.esa.int/archive/ 
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mentioned here, CNEOS 2017-03-09, which have heliocentric 
speeds ~45% and ~19%, respectively, above the hyperbolic 
threshold. 


6. Conclusions 


We presented and analyzed impact data from the meteor 
detected at 2014 January 8 17:05:34 UTC. The U.S. 
Department of Defense verified in a letter referencing this 
work that “the velocity estimate reported to NASA is 
sufficiently accurate to indicate an interstellar trajectory” 
(Shaw 2022). Based on the data provided by CNEOS the 
meteor had an asymptotic speed of v.,~42.1+5.5km s^! 
outside of the solar system. Assuming these parameters 
correspond to the reality, its size, trajectory, and excess speed 
exclude the possibility that it was gravitationally scattered 
within the solar system prior to impact (Wiegert 2014). 

We also analyzed the updated size distribution of interstellar 
meteors, deriving a range of possible slopes of q ~ 3.6 + 0.5, 
consistent with the analytical model by Dohnanyi (1969). 
Finally, we presented a strategy for studying interstellar 
meteors: using a network of ~600 all-sky camera systems to 
determine the orbits and chemical compositions of rz 5 cm 
meteors. This method also allows for the possibility of 
retrieving interstellar meteorite samples. This new field of 
astronomical research is significant as it would save the trip and 
allow us to study samples of materials from other planetary 
systems. 
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